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Abstract. We have analyzed available optical data for Au in the mid-infrared 
range which is important for a precise prediction of the Casimir force. Significant 
variation of the data demonstrates genuine sample dependence of the dielectric 
function. We demonstrate that the Casimir force is largely determined by the 
material properties in the low frequency domain and argue that therefore the 
precise values of the Drude parameters are crucial for an accurate evaluation of 
the force. These parameters can be estimated by two different methods, either by 
fitting real and imaginary parts of the dielectric function at low frequencies, or via 
a Kramers-Kronig analysis based on the imaginary part of the dielectric function 
in the extended frequency range. Both methods lead to very similar results. We 
show that the variation of the Casimir force calculated with the use of different 
optical data can be as large as 5% and at any rate cannot be ignored. To have 
a reliable prediction of the force with a precision of 1%, one has to measure the 
optical properties of metallic films used for the force measurement. 



PACS numbers: 12.20.Ds, 12.20.Fv, 42.50.Lc, 73.61.At, 77.22. Ch 
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1. Introduction 

The Casimir force [1] between uncharged metalhc plates attracts considerable 
attention as a macroscopic manifestation of the quantum vacuum [21 O El HI [6]. 
With the development of microtechnologies, which routinely control the separation 
between bodies smaller than 1 /im, the force became a subject of systematic 
experimental investigation. Modern precision experiments have been performed 
using different techniques such as torsion pendulum [J, atomic force microscope 
(AFM) PIS], microelectromechanical systems (MEMS) [TOl [H HI [H [H [Hj and 
different geometrical configurations: sphere-plate [TUnHHj, plate-plate [Hj and crossed 
cylinders |17| . The relative experimental precision of the most precise of these 
experiments is estimated to be about 0.5% for the recent MEMS measurement [13] 
and 1% for the AFM experiments [9lfT0]. 

In order to come to a valuable comparison between the experiments and the 
theoretical predictions, one has to calculate the force with a precision comparable to 
the experimental accuracy. This is a real challenge to the theory because the force is 
material, surface, geometry and temperature dependent. Here we will only focus on 
the material dependence, which is easy to treat on a level of some percent precision 
but which will turn out difficult to tackle on a high level of precision since different 
uncontrolled factors are involved. 

In its original form, the Casimir force per unit surface [T] 

^^(^) = -24077 W 

was calculated between ideal metals. It depends only on the fundamental constants 
and the distance between the plates L. The force between real materials differs 
significantly from ([T]) for mirror separations smaller than 1 ^m. 

For mirrors of arbitrary material, which can be described by reflection coefficients, 
the force per unit area can be written as [18j : 

d^k f°° dC, [<,k]^e-2KX 



47r2 7o 27r I - r^[iC,,\<.f e-^'^^ 



.= yk2 + J (2) 
where = {rs,rp) denotes the reflection amplitude for a given polarization /i = s, p 



^k2+£«)§+C« 

y^k2+e«)i^-f c^e (*C) 



(3) 



The force between dielectric materials had first been derived by Lifshitz [THl 1^ . 
The material properties enter these formulas via the dielectric function e {iQ at 
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angular imaginary frequencies uj = i^, which is related to the physical quantity 
e" (oj) — Im (e (uj)) with the help of the dispersion relation 

ei^O-l-'-Jd.^. (4) 
n J LU^ + Q'^ 



For metals e" (w) is large at low frequencies, thus the main contribution to the integral 
in Eq. ^ comes from the low frequencies even if Q corresponds to the visible frequency 
range. For this reason the low- frequency behavior of £{lo) is of primary importance. 

The Casimir force is often calculated using the optical data taken from [21j , which 
provides real and imaginary parts of the dielectric function within some frequency 
range, typically between 0.1 and 10** eV for the most commonly used metals, Au, 
Cu and Al, corresponding to a frequency interval [1.519 • 10"'^"', 1.519 • 10^®] rad/s 
(1 eV=1.519 • 10^^ rad/s U). When the two plates are separated by a distance L, 
one may introduce a characteristic imaginary frequency Cch = c/2L of electromagnetic 
field fluctuations in the gap. Fluctuations of frequency C ~ Cch give the dominant 
contribution to the Casimir force. For example, for a plate separation of L = 100 nm 
the characteristic imaginary frequency is Cch = 0.988 eV. Comparison with the 
frequency interval where optical data is available shows that the high frequency data 
exceeds the characteristic frequency by 3 orders of magnitude, which is sufficient for 
the calculation of the Casimir force. However, in the low frequency domain, optical 
data exists only down to frequencies which are one order of magnitude below the 
characteristic frequency, which is not sufficient to evaluate the Casimir force. Therefore 
for frequencies lower than the lowest tabulated frequency, uJc, the data has to be 
extrapolated. This is typically done by a Drude dielectric function 

e{u:) = l~— ^— , (5) 

which is determined by two parameters, the plasma frequency uj-p and the relaxation 
frequency oJt- 

Different procedures to get the Drude parameters have been discussed in the 
literature. They may be estimated, for example, from information in solid state 
physics or extracted form the optical data at the lowest accessible frequencies. The 
exact values of the Drude parameters are very important for the precise evaluation of 
the force. Lambrecht and Reynaud [18_ fixed the plasma frequency using the relation 

where N is the number of conduction electrons per unit volume, e is the charge and 
m* is the effective mass of electron. The plasma frequency was evaluated using the 
bulk density of Au, assuming that each atom gives one conduction electron and that 
the effective mass coincides with the mass of the free electron. The optical data at the 
lowest frequencies were then used to estimate u>r with the help of Eq. ([5]). In this way 
the plasma frequency Wp — 9.0 eV and the relaxation frequency = 0.035 eV have 
been found. This procedure was largely adopted in the following [SJ [TTl HUl [IHl E] ■ 
However, on the example of Cu, it was stressed in |18j that the optical data may vary 

X In |18) a conversion factor 1.537 • 10^^ rad/s was used, leading however to a negligible difference in 
the Casimir force (well below 1%). 
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from one reference to another and a different choice of parameters for the extrapolation 
procedure to low frequencies can influence the Casimir force significantly. 

Bostrom and Sernelius P3' and Svetovoy and Lokhanin extracted the low- 
frequency optical data by fitting them with Eq. ([5]) . For one set of data from Ref . [53] 
the result [22] was close to that found by the first approach, but using different sources 
for the optical data collected in Ref. [25] an appreciable difference was found [24l [22] . 
This difference was attributed to the defects in the metallic films which appear as the 
result of the deposition process. It was indicated that the density of the deposited 
films is typically smaller and the resistivity larger than the corresponding values for 
the bulk material. The dependence of optical properties of Au films on the details 
of the deposition process, annealing, voids in the films, and grain size was already 
discussed in the literature |26] . 

In this paper we analyze the optical data for Au from several available sources, 
where the mid-infrared frequency range was investigated. The purpose is to establish 
the variation range of the Drude parameters and calculate the uncertainty of the 
Casimir force due to the variation of existing optical data. This uncertainty is of great 
importance in view of the recent precise Casimir force measurement |27l[T3] which have 
been performed with high experimental accuracy. On the other hand, sophisticated 
theoretical calculations predict the Casimir force at the level of 1% or better. These 
results illustrate the considerable progress achieved in the field in only one decade. 
In order to assure a comparison between theory and experiment at the same level of 
precision, one has to make sure that the theoretical calculation considers precisely the 
same system investigated in the experiment. This is the key point we want to address 
in our paper. With our current investigation we find an intrinsic force uncertainty 
of the order of 5% coming from the fact that the Drude parameters are not precisely 
known. These parameters may vary from one sample to another, depending on many 
details of the preparation conditions. In order to assure a comparison at the level of 
1% or better between theoretical predictions and experimental results for the Casimir 
force, the optical properties of the mirrors have to be measured in the experiment. 

The paper is organized as follows. In Sec. [5]we explain and discuss the importance 
of the precise values of the Drude parameters. In Sec. [3] the existing optical data for 
gold are reviewed and analyzed. The Drude parameters are extracted from the data 
by fitting both real and imaginary parts of the dielectric function at low frequencies in 
Sec. m In Section [5] the Drude parameters are estimated by a different method using 
Kramers-Kroning analysis. The uncertainty in the Casimir force due to the sample 
dependence is evaluated in Sec. [S]and we present our conclusions in Sec. [7| 

2. Importance of the values of the Drude parameters 

In Figure [T] (left) we present a typical plot of the imaginary part of the dielectric 
function, which comprises Palik's Handbook data for gold jH]. The solid line shows 
the actual data taken from two original sources: the points to the right of the arrow 
are those by Theye [28] and to the left by Dold and Mecke [29]. No data is available 
for frequencies smaller than the cutoff frequency uJc (0.125 eV for this data set) and 
e" (w) has to be extrapolated into the region oj < ujc- The dotted line shows the 
Drude extrapolation with the parameters tOp = 9.0 cV and lOt- = 0.035 eV obtained in 
Ref. [IH]. 

One can separate three frequency regions in Fig. [T] (left panel). The region marked 
as 1 corresponds to the frequencies smaller than lOc- The region 2 defining the Drude 
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Figure 1. Left panel: Palik's Handbook data for Au 21 (solid line) extrapolated 
to low frequencies (dotted line) with the Drude parameters indicated in the corner. 
Right panel: contributions of different real frequency domains to the dielectric 
function on the imaginary axis e{iQ. 



parameters extends from the cutoff frequency to the edge of the interband absorption 
ujo- The high energy domain ui > ujq is denoted by 3. 

We may now deduce the dielectric function at imaginary frequencies ^ using the 
Kramers-Kronig relation 

e(*C) = l+£i(*C)+e2(«C)+e3(iC), (7) 

where the indices 1,2, and 3 indicate respectively the integration ranges < lu < lOc, 
LOc < oj < ojQ, and wq < w < oo. ei can be derived using the Drude model ([5]) leading 
to 

The two other functions £2 and £3 have to be calculated numerically. The results 
for all three functions as well as for £ (iC) are shown in Fig. [1] (right). One can 
clearly see that £1 dominates the dielectric function at imaginary frequencies up 
to C ~ 5 eV. £2 (iC) gives a perceptible contribution to e , while £3 (iC) produces 
minor contribution negligible for C < 0.5 eV. 

As mentioned in the Introduction, we may introduce a characteristic imaginary 
frequency Cch = c/2L of field fiuctuations which give the dominant contribution 
to the Casimir force between two plates at a distance L. For a plate separation 
oi L = 100 nm the characteristic imaginary frequency is — 0.988 eV. At this 
frequency the contributions of different frequency domains to e (iCch) are £1 — 68.42, 
£2 — 15.65, and £3 — 5.45. This means that for all experimentally investigated 
situations, L > 100 nm, region 1, corresponding to the extrapolated optical data, 
gives the main contribution to s It is therefore important to know precisely the 
Drude parameters. 
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Figure 2. Left panel: Available optical data in the mid-infrared region. The 
dots represent the Dold and Mecke data for lj < I eV I29| and Theye data |28j for 
higher frequencies. The squares denote the Weaver data 1301 . The circles stand 
for the data from |40) . The triangles represent the data [41]. Solid squares, circles, 
and triangles are used to mark e" (ut) while the open symbols are used for e' (lu). 
Right panel: e"{Lj) in the interband region for different samples. The solid line 
represents the data measured with the well annealed bulk-like film by Theye 1281 . 
The dots are the data by Johnson and Christy |37| found for unannealed films. 
The dashed and dash-dotted lines are recent data sets by Wang et al. |34| for 
unannealed films. They correspond to films deposited with e-beam and thermal 
evaporation methods, respectively. 



3. Analysis of different optical data for gold 

The optical properties of gold were extensively investigated in 50-70th. In many 
of those works the importance of sample preparation methods was recognized and 
carefully discussed. A complete bibliography of the publications up to 1981 can be 
found in Ref . [30j . Regrettably the contemporary studies of gold nanoclusters produce 
data inappropriate for our purposes. Among recent experiments let us mention the 
measurement of normal reflectance for evaporated gold films [31) , which was performed 
in the wide wavelength range 0.3 — 50 fim, but unfortunately does not permit to 
evaluate independently both real and imaginary parts of the dielectric function. In 
contrast, the use of new ellipsometric techniques [5^ [55] has produced data for the 
real and imaginary part of the dielectric function for energy intervals 1.5 — 4.5 eV [33] 
and 1.5-3.5 eV [35 . 

A significant amount of data in the interband absorption region (domain 3) has 
been obtained by different methods under different conditions [36l [28l jST] |38l (39] 
[Ml [55] . Though this frequency band is not very important for the Casimir force, it 
provides information on how the data may vary from one sample to another. On the 
contrary there are only a few sources where optical data was collected in the mid- 
infrared (domain 2) and from which the dielectric function can be extracted. The 
data available for e' (w) and e" (w) in the range lu < 1.5 eV and interband absorption 
domain 3 are presented respectively in the left and right graph of Fig. |2j These data 
sets demonstrate considerable variations of the dielectric function from one sample to 
another. 
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Let us briefly discuss the sets of data [13 [301 HOI HI] used in our analysis and 
the corresponding samples. The commonly used Handbook of Optical Constants of 
Solids [H] comprises the optical data covering the region from 0.125 to 9184 eV (dots 
in Fig. [2). The experimental points are assembled from several sources. For w < 1 eV 
they are reported by Dold and Mecke [29]. For higher frequencies up to 6 eV they 
correspond to the Theye data [28] . Dold and Mecke give only little information about 
the sample preparation, reporting that the films were evaporated onto a polished glass 
substrate and measured in air by using an ellipsometric technique 29J. Annealing of 
the samples was not reported. 

Theye |28| described her films very carefully. The samples were semitransparent 
Au films with a thickness of 100 — 250 A evaporated in ultrahigh vacuum on 
supersmooth fused silica. The substrate was kept in most cases at room temperature. 
After the deposition the films were annealed in the same vacuum at 100 — 150° C. The 
structure of the films was investigated by X-ray and transmission-electron-microscopy 
methods. The dc resistivity of the films was found to be very sensitive to the 
preparation conditions. The errors in the optical characteristics of the films were 
estimated on the level of a few percents. 

The handbook [30] embraces the optical data from 0.1 eV to 28.6 eV (marked 
with squares in Fig. 2). The data in the domain oj < 4 eV is provided by Weaver et 
al. [3D]. The values of e{uj) were found for the electropolished bulk Au(llO) sample. 
Originally the refiectance was measured in a broad interval 0.1 < w < 30 eV and then 
the dielectric function was determined by a Kramers-Kronig analysis. Due to indirect 
determination of e the recommended accuracy of these data sets is only 10%. 

The optical data of Motulevich and Shubin [40] for Au films is marked with circles 
in Fig. 2. In this paper the films were carefully described. Gold was evaporated on 
polished glass at a pressure of ^ 10~^ Torr. The investigated films were 0.5 — 1 /im 
thick. The samples were annealed in the same vacuum at 400° C for more than 3 
hours. The optical constants n and k {n + ik = y^) were measured by polarization 
methods in the spectral range 1 — 12 fim. The errors in n and k were estimated as 
2-3% and 0.5-1%, respectively. 

Finally, the triangles represent Padalka and Shklarevskii data [|41 for unannealed 
Au films evaporated onto glass. 

The variation of the data points from different sources cannot be explained by 
experimental errors. The observed deviation is the result of different preparation 
procedures and refiects genuine difference between samples. The deposition method, 
type of the substrate, its temperature, quality and the deposition rate infiuence the 
optical properties. When we are speaking about a precise comparison between theory 
and experiment for the Casimir force at the level of 1% or better, there is no such 
material as gold in general any more. There is only a gold sample prepared under 
definite conditions. 

4. Evaluation of the Drude parameters through extrapolation 

We will now use the available data in the mid-infrared region to extrapolate into the 
low frequency range. If the transition between inter- and intraband absorption in gold 
is sharp, the data below ujq should be well described by the Drude function 



Sample dependence of the Casimir force 



8 



For Lu ^ lUt, the data on the log- log plot should fit straight lines with the slopes 
—2 and —3 for e' and e", respectively, shifted along the ordinate due to variation of 
the parameters for different samples. The data points in the right graph of Fig. [5] 
are in general agreement with these expectations. The onset values for e", hi^ujpUJr), 
vary more significantly due to a significant change in LUr for different samples, but the 
Casimir force is in general not very sensitive to the relaxation parameter jl8l. The 
onset values for — e', In(wp), vary less but this variation is more important for the 
Caimir force, which is particularly sensitive to the value of the plasma frequency uip. 
The Drude parameters can be found by fitting both e' and e" with the functions 
This procedure is discussed below. 

The dielectric function for low frequencies, lu < uJc, is found by the extrapolation 
of the optical data from the mid- infrared domain, tOc < uj < ujq- The real and 
imaginary parts of e follow from Eq. ([9]) with an additional polarization term V 
in e': 

^'H = ^-^^' ^"i-)- (10) 

The polarization term appears here due to the following reason. The total dielectric 
function e — £(c) -t- includes contributions due to conduction electrons e(c) and the 
interband transitions £(i). The polarization term consists of the atomic polarizability 
and polarization due to the interband transitions e'^^^ 

p = l + ^+,' (^), (11) 
£o ^ ' 

where a is the atomic polarizability and Na the concentration of atoms. If the 
transition from intra- to interband absorption is sharp, the polarization can be 
considered as constant, because the interband transitions have a threshold behavior 
with an onset frequency loq and the Kramers-Kronig relation allows one to express e'^^j 
as 

^«H— y^-^- (12) 

For w <C Wo this integral does not depend on w, leading to a constant e'^^^ (w). In 
reality the situation is more complicated because the transition is not sharp and many 
factors can influence the transition region. We will assume here that 7-" is a constant 
but the fitting procedure will be shifted to frequencies where the transition tail is not 
very important. In practice Eq. (jlOp can be applied for ld < 1 eV. 

Our purpose is now to establish the magnitude of the force change due to 
reasonable variation of the optical properties. To this end the available low-frequency 
data for e' [lu) and e" [lo) presented in the left graph of Fig. [5] were fitted with Eq. 
pop. The results together with the expected errors are collected in Table [TJ 

The error in Table[T]is the statistical uncertainty. It was found using a criterion 
for joint estimation of 3 parameters [33j- For a given parameter the error corresponds 
to the change A^^ = 1 when two other parameters are kept constant. The parameter 
V enters (|10|) as an additive constant and in the considered frequency range its value 
is smaller than 1% of e' {lo) . That is why the present fitting procedure cannot resolve 
it with reasonable errors. 
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N 


ujp{eV) 


OJr ■ 102 (eV) 


V 




1 


7.50 ±0.02 


6.1 ±0.07 


-27.67 ±5.79 


Palik, 66 points , • 


2 


8.41 ±0.002 


2.0 ±0.005 


7.15 ±0.035 


Weaver, 20 points, ■, □ 


3 


8.84 ±0.03 


4.2 ±0.06 


12.94 ± 16.81 


Motulevich, 11 points, 9,0 


4 


6.85 ±0.02 


3.6 ±0.05 


-12.33±9.13 


Padalka 11 points, T, V 



Table 1. The Drude parameters found by fitting the available infrared data for 
e' {ui) and e" (ui) with Eq. mOI I. The error is statistical. 



As mentioned before, in the case of the Weaver data [3D] the recommended 
precision in e' and e" is 10% while Motulevich and Schubin reported 2-3% and 0.5-1% 
errors in n and k. We did not take these errors exphcitly into account as we do not 
know if they are of statistical or systematic nature or a combination of both. But 
to illustrate their possible influence let us just mention that if we interpret them as 
systematic errors, we can propagate the errors in e or 71, k to the values of LUp and uJr, 
leading to an additional error in Up of about 5% for the Weaver data and 1% for the 
Motulevich data and twice as large in ujr- 

Significant variation of the plasma frequency, well above the errors, is a distinctive 
feature of the table. The bulk and annealed samples (rows 2 and 3) demonstrate larger 
values of Wp. The rows 1 and 4 corresponding to the evaporated unannealed films give 
rise to considerably smaller plasma frequencies Up. Note that our calculations are in 
agreement with the one given by the authors [521 HI] themselves. 

To have an idea of the quality of the fitting procedure, we show in Fig. [3| the 
experimental points and the best fitting curves for Dold and Mecke data [29l [21| (full 
circles and solid lines) and Motulevich and Shubin data HO] (open circles and dashed 
lines). Only 25% of the points from [5T] are shown for clarity. One can see that for 
e" at high frequencies the dots lie above the solid line demonstrating presence of a 
wide transition between inter- and intraband absorption. Coincidence of the solid and 
dashed lines for e" is accidental. The fits for e' are nearly perfect for both data sets. 

It is interesting to see on the same figure how well the parameters u!p = 9.0 eV, 
lOt = 0.035 eV agree with the data in the mid-infrared range. The curves corresponding 
to this set of parameters are shown in Fig. [3] as dotted lines. One can see that the 
dotted line, which describes e" is very close to the solid line. However, the dotted 
line for e' does not describe well the handbook data (full circles). It agrees much 
better with Motulevich and Shubin data [40] (open circles). The reason for this is 
that LOp = 9.0 eV is the maximal plasma frequency for Au. Any real film may contain 
voids leading to smaller density of electrons and, therefore, to smaller ojp. Motulevich 
and Shubin [JD] annealed their films which reduced the number of defects and made 
the plasma frequency close to its maximum. A plasma frequency Wp = 9.0 eV was 
also reported in Ref. [44] , where the authors checked the validity of the Drude theory 
by measuring refiectivity of carefully prepared gold films in ultrahigh vacuum in the 
spectral range 0.04 < a; < 0.6 eV. Therefore, this value is good if one disposes of well 
prepared samples. 

5. The Drude parameters from Kramers-Kronig analysis 

Because the values of the Drude parameters are crucial for a reliable prediction of 
the Casimir force, it is important to assess that different methods to determine the 
parameters give the same results. Alternatively to the extrapolation procedure of the 
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Figure 3. The infrared optical data by Dold and Mecke |29| (full circles) and 
by Motulevich and Shubin 40, (open circles) together with the best Drude fits 
given by the solid and dashed lines, respectively. The dotted lines present the fit 
with LUp = 9 eV and lj^- = 35 meV which agrees better with the Motulevich and 
Shubin data (open circles) than with the handbook data (full circles). 



previous section we will now discuss a procedure based on a Kramers-Kronig analysis. 
To this aim we will extrapolate only the imaginary part of the dielectric function to 
low frequencies uj < Uc- The dispersion relation between e' and e" 



can then be used to predict the behavior of e'{uj) and compare it with the one observed 
in the experiments. From this comparison the Drude parameters can be extracted. 

The low- frequency behavior of s"{lu) is important for the prediction of e' because 
for metals e"{Lu) 3> 1 in the low frequency range. Therefore, at w < we are 
using £"(0-1) from Eq. At higher frequencies the experimental data from different 
sources [HI |30l |40l [41] are used. The data in Refs. [40l [41] must be extended to high 
frequencies starting from cu = 1.25 eV. We do this using the handbook data [21]. 

Let us start from the data for bulk Au(llO) ^0]. This data set is given in 
the interval 0.1 < w < 30 eV. Below w = 0.1 eV we use the Drude model for 
e" and above uj = 30 eV the cubic extrapolation C/uj^. The Drude parameters 
are practically insensitive to the high frequency extrapolation. The data set was 
divided into overlapping segments containing 12 points. Each segment was fitted with 
a polynomial of forth order in frequency. The first segment, were e"(w) increases very 
fast, was fitted with the polynomial in l/uj. Then, in the range of overlap (4 points) a 
new polynomial smoothly connecting two segments was chosen. In this way we have 
fitted the experimental data with a function which is smooth up to the first derivative. 

The real part of the dielectric function e'{uj) is predicted by Eq. ([T3|) as a function 
of the Drude parameters tOp and LOr ■ These parameters are chosen such as to minimize 
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Figure 4. |e'| as a function of for bulk gold. Dots are the experimental 
data I30| . The solid line is the prediction according to Eq. (I13I I with the Drude 
parameters ijp = 8.40 eV, oj^ = 0.02 eV . 



the difference between observed and predicted values of e'(tj), leading to Wp — 8.40 eV 
and u>T — 0.020 eV. These parameters are in reasonable agreement with the ones 
indicated in Tab. [TJ In Fig. [4] the experimental data (dots) and |e'('j-')| found from 
Eq. (jl3p (solid line) are plotted, showing perfect agreement at low frequencies, while 
at high frequencies lo > 2.6 eV the agreement is not very good. This may be fixed by 
choosing an appropriate high frequency extrapolation. We do not give these details 
here as this extrapolation has practically no influence on the Drude parameters. 

When applying the same procedure to the handbook data [21], we find ujp = 
7.54 eV and lOt = 0.051 eV, again in agreement with the parameters indicated in 
Tab.[T] Fig.[5]shows a plot of e'{uj) predicted with these parameters. At low frequencies 
the agreement with the experimental data is good but it becomes worse when the 
interband data joins the intraband (high frequency) data [28j. These two data 
sets correspond to samples with different optical properties. In this case the dispersion 
relation (fT3|) is not necessarily very well satified. In contrast with the previous case, 
high frequency extrapolation cannot improve the situation; it influences the curve only 
marginally. 

Following the same procedure for the Motulevich and Shubin data [40], we find 
the Drude parameters cUp = 8.81 eV, uJt = 0.044 eV which are close to the values in 
Tab. [TJ The experimental data and calculated function |e'(w)| are shown in Fig. [Hi 
There is good agreement for frequencies a; < 4 eV as the data in Ref. [ID] matches very 
well the Theye data [15] • Deviations at higher frequencies are again quite sensitive to 
high-frequency extrapolation as already noted before. 

Similar calculations done for the Padalka and Shklyarevskii data [41] give the 
Drude parameters tUp = 6.88 eV and uJt = 0.033 eV, producing good agreement only 
in the range lu < 1.3 eV because this data set matches only poorly the Theye data 





Figure 6. |e'| as a function of oj for Motulevich and Shubin data |40) extended 
by the handbook data 21 for uj > 1.25 eV (dots). The sohd line is found from 
Kramers-Kronig relation. 
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Using the Kramers-Kronig analysis for the determination of the Drude parameters 
leads essentially to the same parameters for all 4 sets of the experimental data. 
Experimental and calculated curves for s'{uj) are in very good agreement at low 
frequencies. At high frequencies the agreement is not so good for two different 
reasons. First, at high frequencies the calculated curve is sensitive to the high- 
frequency extrapolation and thus a better choice of this extrapolation can significantly 
reduce high frequency deviations. The other reason is that one has to combine the 
data from different sources to make a Kramers-Kronig analysis possible. These data 
sets do not always match each other well as it is for example the case of the Dold and 
Mecke data and the Theye data. In this case significant errors might be introduced 
in the dispersion relation. Indeed the Kramers-Kronig analysis is a valuable tool only 
for data taken from the same sample. 



6. Uncertainty in the Casimir force due to variation of optical properties 

We will now assess how the values of the Casimir force arc influenced by the different 
values of the Drude parameters. As an example we consider as input the optical data 
for Au from [2T. 

Instead of calculating the absolute value of the Casimir force, we will give the 
factor which measures the reduction of the Casimir force with respect to the ideal 
Casimir force between perfect mirrors as introduced in |18| 



J dnn'J dCY.-^^^, (14) 



VF 



The dielectric function at imaginary frequencies e{iQ is calculated using the Kramers- 
Kronig relation ([¥]) and the integration region is divided in two parts 

xe"{x) , r r--) xe"{x) , , r 

We assume that for x < Xc the Drude model ^ is applicable. Then the integration 
in Ii may be carried out explicitly, see ([H]). In I2 we integrate from Xc = 0.125 eV to 
Xma,x = 9000 eV (corresponding to the range of available optical data in [21]). 

For the calculation of the reduction factor p4|) the integration range was chosen 
as 10""* — 10'^ eV. We also varied the integration range by half an order of magnitude, 
which changed the result by less than 0.1%. The results of the numerical integration 
are collected in Table [2J 

The first four rows of the table present the reduction factors for four pairs of 
the Drude parameters that were obtained by fitting the optical data from different 
sources. The next row shows the result obtained for cjp = 9 eV and oJt — 35 meV. 
The last two rows show the variation of the reduction factor if the plasma frequency 
ujp or the relaxation parameter uJt are varied by ±15% and ±30%, respectively. The 
upper (lower) line corresponds here to the upper (lower) sign. 

The variation of the optical data and the associated Drude parameters introduces 
a variation in the Casimir force ranging from 5.5% at short distances (100 nm) to 1.5% 
at long distances (3 ^m). The distance dependence is of course related to the fact 
that the material properties influence the Casimir force much more at short than 
at long plate separation. The strongest variation of 5.5% gives an indication of the 
genuine sample dependence of the Casimir force. For this reason it is necessary to 
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Up, ujr{eV) \L(^m) 


0.1 


0.3 


0.5 


1.0 


3.0 


1. 


ujp = 7.50, Ur = 0.061 


0.43 


0.66 


0.75 


0.85 


0.93 


2. 


uj^ = 8.41, ujr = 0.02 


0.45 


0.69 


0.79 


0.88 


0.95 


3. 


= 8.84, = 0.0422 


0.46 


0.69 


0.78 


0.87 


0.94 


4. 


Ljp — 6.85, = 0.0357 


0.42 


0.65 


0.75 


0.84 


0.93 


5. 


= 9.00, ujr = 0.035 


0.47 


0.71 


0.79 


0.88 


0.95 


6. 


t^p = 7.50 ± 15% 
UJr = 0.061 


0.45 
0.41 


0.68 
0.63 


0.77 
0.73 


0.86 
0.83 


0.94 
0.92 


7. 


wp = 7.50 

= 0.061 ± 30% 


0.42 
0.44 


0.65 
0.67 


0.74 
0.76 


0.84 
0.86 


0.92 
0.93 



Table 2. The reduction factors at different plate separations calculated with 
the different pairs of values of the Drude parameters corresponding to different 
data. The last two rows show the variation of the reduction factor when either 
the plasma frequency or the relaxation parameter is varied. 

measure the optical properties of the plates used in the Casimir force measurement if 
a precision of the order 1% or better in the comparison between experimental values 
and theoretical predictions is aimed at. Incidentally let us notice that the plasma 
frequency uj-p — 7.5 eV, which is found here to fit best Palik's handbook data |21j . 
is basically the same as the one proposed alternatively in [T5] for Cu, which has 
very similar optical properties to Au concerning the Casimir force [45]. For Cu, the 
variation of the plasma frequency from Wp = 9 eV to ujp = 7.5 eV introduced a 
variation of the Casimir force up to 5% |Tg] . 

In order to asses more quantitatively the role of the two Drude parameters, we 
show in the last two rows of table [5] the variation of the reduction factor when either 
the plasma frequency or the relaxation parameter is varied with the other parameter 
kept constant. One can see that the increase (decrease) of the relaxation parameter 
by StOr = 30% lowers (increases) the reduction factor rjp &t L = 0.1 iim by only 
SrjF = 1.6%. However, the 15% variation of the plasma frequency leads to 4.2% 
change in the reduction factor. Thus the Casimir force is much more sensitive to the 
variation of the plasma frequency, basically as the plasma frequency determines the 
reflection quality of the plates (an infinite plasma frequency corresponds to perfectly 
reflecting mirrors). 

7. Conclusions 

In this paper we have performed the first systematic and detailed analysis of optical 
data for Casmir force measurements. We have studied the relative importance of the 
different frequency regions for the Casimir force as a function of the plate separation 
and established the critical role of the Drude parameters in particular for short 
distance measurements. We have then analyzed and compared four different sets 
of optical data. For each set we have extracted the corresponding plasma frequency 
and relaxation parameter either by fitting real and imaginary part of the dielectric 
function at low frequencies or by using a detailed Kramers-Kronig analysis. Both 
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methods lead essentially to the same results. The Kramers-Kronig analysis reveals 
itself to be a powerful tool for the estimation of the low frequency Drude parameters 
for data coming from the same sample. 

A variation of the values of the Casimir force up to 5.5% is found for different 
optical data sets. This gives an intrinsic unknown parameter for the Casimir force 
calculations and demonstrates the genuine sample dependence of the Casimir force. 
The today existing numerical and analytical calculations of the Casimir force in 
themselves are very precise. In the same way, measurements of the Casimir force 
have achieved high accuracy over the last decade. In order to compare the results 
of the achievements in theory and experiment at a level of 1% precision or better, 
the crucial point is to make sure that calculations and experiments are performed 
for the same physical sample. One therefore has to know the optical and material 
properties of the sample used in the experiment. These properties must be measured 
for frequencies as low as possible. In practice, the material properties have to be known 
over an interval of about 4 orders of magnitude around the characteristic frequency 
Cch = c/2L. For a plate separation of L = 100 nm this means an interval [10 meV, 
100 eV]. If measurements at low frequencies arc not possible, the low frequency Drude 
parameters should be extracted from the measured data, by one of the two methods 
discussed here. 
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